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ABSTRACT 

Prom common proper motion and signatures of youth, researchers have identified 
about 30 members of a putative TW Hydrae Association. Only four of these had 
parallactic distances from Hipparcos. We have measured parallaxes and proper motions 
for 14 primary members. We combine these with literature values of radial velocities 
to show that the Galactic space motions of the stars, with the exception of TWA 9 and 
22, are parallel and do not indicate convergence at a common formation point sometime 
in the last few million years. The space motions of TWA 9 and 22 do not agree with 
the others and indicate that they are not TWA members. The median parallax is 18 
mas or 56 pc. We further analyze the stars' absolute magnitudes on pre-main sequence 
evolutionary tracks and find a range of ages with a median of 10.1 Myr and no correlation 
between age and Galactic location. The TWA stars may have formed from an extended 
and filamentary molecular cloud but are not necessarily precisely coeval. 

Subject headings: open clusters and associations: individual (TW Hydrae) - stars: 
distances - stars: kinematics and dynamics - stars: pre-main sequence 



1. Introduction 

Nearby young stars of age ~5-10 Myr provide our best opportunity to study the late stages of 
star and planet formation with high sensitivity and spatial resolution. During this time period, the 
last gas-rich disks dissipate, and the onset of the debris disk phase oc curs. The star TW Hydrae 



sports a massive disk and was first identified as an isolated T Tauri star (jRucinski &: Krautterll 19831 ) 
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and then the first member of an association of young stars (|de la Reza et al.lll989l : lGregorio-Hetem et al 
1992; iKastner et alj 119971 ) . As the nearest example of a protoplanetary disk, TW Hya has been 



studied at every available wavelength and spatial resolution, but to understand its disk in context, 
TW Hya's age must be well-determined and its disk evolution compared to other stars of similar 
age and mass. 

Over the last decade, ~30 members of a putative TW Hydrae Association (TW A) have been 



ident i fied from a combination of stellar-activity searches and s pace motion stud i es fe.g. IWebb et al 



1999 



Sterzik et al 



1999 



Zuckerman et alJ l200ll : iGizisI l2002l : ISong et alJ l2003l : IScholz et al 



2005 



Looper et al.ll2007l : IShkolnik et al.ll201ll ). Members range in spectral type from AO to brown dwarfs 
and have presumed common ages ~10 Myr, as determined from L i depletion and modeling on 
pre-main sequence evolutionary tracks (jBarrado Y Navascued l2006l ) . Disks in TWA range from 
the four accreting, gas rich, protoplanetary ones (TWA 1, 3, 27, and 30), to seven transitional 
or debris disks (TWA 4B, 7, 11 A, 26, 28, 31 and 32), to the majority of inembers that have no 



2009 



detectable disks at all (|Weinberger et alJl2004l : ILow et al. 



2005 



Riaz &: Gizis 



2008 



Plavchan et al 



Schneider et al. 2OI2I ). 



The Hipparcos mission determined distances to only four TWA members, including TW Hya 
itself, for an average distance of 70 pc. Although the association contains the youngest stars close 
to the Sun, most of its members are low-mass stars and were therefore too faint for Hipparcos. 
Ground-b ased parallaxes ha ve been measured for two bro wn dwarfs in the as sociation - TWA 27 



(2M1207) (iGizis et al.ll2007l ^ and TWA 28 (SSSPM ,71102) (ITeixeira et a 
star that may not actually be a member - TWA 22AB, (jTeixeira et al. 



20081 ) , and one additional 



20091). This highlights the 



problem of determining association membership - although their proper motions are similar, as 
ensured by search methods, without distances and therefore true Galactic space motions TWA 
stars can only be presumed to be associated. Hence, various members have been suggested and 
then excluded by proper rnotion studies using varyingly restr ictive conditions for membership (e.g. 



Makarov &: Fabricius 



2001 



Song et al.ll2003l : lMamaiekll2005l ). 



Distances to a substantial number of members would aid in defining the association and es- 
timating its age from pre-main sequence tracks. If the kinematics allow, a cluster expansion age 
could also be determined. In this paper we present parallaxes and proper motions to 14 primary 
stars and two visual binary companions. 



Data 



We have observed 14 TWA primary members with the CAPSCam instrument at the 2.5m 
DuPont Telescope at L as Campan a s Obs ervat ory. The instrument and bas ic data reduction tech- 
niques are described in Iboss et al.l (2009) and Anglada-Escude et al. ( 2012 ) and only briefly sum- 
marized here. CAPSCam utilizes a Hawaii-2RG HyViSI detector filtered to a bandpass of 100 nm 
centered at 865 nm with 2048x2048 pixels, each subtending 0'.'196 on a side. The main advantage 
of this camera is its ability to achieve simultaneously high S/N on a bright target star and fainter 
astrometric reference stars. Our typical TWA sources have I~10 (see Table[T|); we place these bright 
target stars in an independently readable subarray called the "guide window" (GW). We typically 
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locate the GW in the center of the full field (FF) and rapidly read and reset it for integration times 
down to 0.2 s for a 64x64 pixel subarray. Simultaneously, the FF integrates for as long as necessary 
to obtain high S/N on many reference stars. Table [1] gives typical integration times for each target 
in the GW and FF, although the times were adjusted during each epoch to account for seeing and 
clouds. 

We drew our sample from stars without parallaxes (thus excluding TWA 1, 4AB, 9AB, llA, 
22, 27, aiid 28). We only observed bona fide members identified in the coi ivergent point an alysis of 



Mamaiekl ((20051); thus we excluded TWA 17, 18, 19AB, and 24. Although iMamaiekl (120051) ex clude 
TWA 12, we included it because of its large discrepancy between photometric (jSong et al.ll2003l ) 
and kinematic distance. We excluded TWA 3AB for being too bright. CAPSCam saturates in 
0.2 s on I~9. We excluded TWA 6, 7, and IIB for being near very bright stars that would fall 
within the CAPSCam field of view because saturated images leave long-lasting after-images on the 
detector. We excluded TWA 8AB because the two stars, separated by 13" do not fit in the 64x64 



GW. Three targets - 30 AB, 31, and 32, were discovered after our survey began (ILooper et al 



irgets - 

Shkolnik et aPboill ). The remaining stars were all observed with the exception of TWA 10 (Table 



201C 



[T]). CAPSCam could in the future provide parallaxes for the remaining stars with I>9 - TWA 8AB, 
10, 17, 18, 30, 31 and 32. 

We can operate in a mode, called "Guide-window shutter" (GWS) where an iris shutter opens 
over the FF only when the GW is integrating and closes during the initial pixel resets and during 
each GW read. This is done to guarantee that the images in the GW and FF sample the atmosphere 
identically and that no astrometric bias is introduced between them. Some of our data were taken 
with the shutter in this mode, but some were taken with no shutter (NS) and some in a mode where 
the shutter closes during the initial pixel resets but remains open during each GW read (NGWS). 
In principle, the best astrometric precision is obtained in GWS mode, but we operated without it 
in order to improve efficiency. 

Flat field images to correct for pixel-to-pixel variations either intrinsic to the detector or due 
to the finite opening time of the iris shutter are obtained while exposing on a quartz lamp projected 
on a flat-field screen. 

The imaging of each astrometric fleld is repeated typically 15-20 times at each epoch. For each 
FF image, we obtain many GW images. In post-processing, these are summed and and re-inserted 
into the FF image. For each night, one image is selected as an astrometric template of the fleld for 
that night. Sources are found automatically, and a flne centroiding algorithm is applied to produce 
a catalog with the sub-pixel positions of all of the objects in all of the images of the fleld for a given 
night. 

Each star was observed in at least four independent epochs, with all but one star observed in 
flve or more, i.e. dates separated by enough days that they provide independent constraints on 
the parallax. The dates of observation are given in Tabled! Data from all epochs are combined in 
an astrometric solution to derive the positions, proper motions, and parallaxes of all the stars in 
each target fleld. The astrometric solution is an iterative process. An initial catalog of positions is 
generated from a chosen epoch, a transformation is applied to every other epoch's catalog to match 
the initial catalog, and the apparent trajectory of each star is then fltted to a basic astrometric 
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model. The initial catalog is updated with new positions, proper motions, and parallaxes, and a 
subset of well-behaved stars is selected to be used as the reference frame. The reference stars must 
be successfully extracted in every epoch and a subset of at least 15, and more typically 30, are 
chosen that show the smallest epoch-to-epoch variation in their solutions. Over all of our target 
fields, these stars have typical I-band magnitudes of 13.5 - 17.9 with a media n of 16.2. This process 



i s the n iterated a small number of times. Again, details may be found in lAnglada-Escude et al 



2.1. Zero-Point Parallax Correction 

The motion of the target star is measured with respect to background stars, which are not 
truly stationary and that all have parallactic motions given by the Earth's motion and therefore 
move in the same direction. This introduces a small bias, so the average parallax of the reference 
stars must be removed to find the absolute parallax. 

If all the stars in the field had perfectly known distances and that information was inserted a 
priori, the parallax zero-point correction would always be a positive number. Note, however, that 
some zero-point corrections in Table [2] are negative. Because we do not know the distances a priori, 
in the astrometric solution, the parallaxes for all objects are initialized to zero. In the iterations 
that follow, each individual parallax and proper motion is adjusted. While the mean parallax 
measurement over all the stars after the first iteration should still be approximately zero, the mean 
parallax of a subset of them, i.e. those used as reference stars, can be either positive or negative 
due to statistical fluctuations. At any epoch, the position of a star has centroiding uncertainties, 
and for distant stars, proper motion will take out all apparent motion of the star leaving positional 
residuals that are both positive and negative. Therefore, although the true parallax to every star 
must be positive, we allow the fit parallaxes to take on positive and negative values. The quality 
of the final astrometric solution as measured by the residuals on all the stars is independent of the 
value of the mean zero-point parallax and will not be adjusted in subsequent iterations. 

To find the zero-point for each field, we estimate a photometric distance to the br ightest refer- 



ence stars by fitting a Kurucz stella r model to catalogued USNO-Bl - B2, R2, and I; (IMonet et al 



20031 ) and 2MASS - J, H, and (jSkrutskie et al.ll2006l ) photometry and assuming each star is a 



dwarf. Giant stars are then easily recognizable because they appear to be so close that they should 
have detectable parallaxes, and we refit them as giants. Dwarf stars with fit Toff <3800 K are not 
considered because the stellar models are less reliable. We average the difference between our astro- 
metrically determined (even if they are not statistically significant) and photometric parallaxes to 
find the average bias and its uncertainty and subtract it from our relative parallaxes and propagate 
th e uncertainty. A comparison of our parallaxes determined this way to literature values is given 



m 



Anglada-Escude et al.l ((2012) 



In principle, this zero-point correction could be done for the proper motions as well, but 
the reference stars do not generally have catalogued proper motions to use in measuring their bias. 
Instead, we estimate our bias directly by comparing the proper motions of th e TWA stars themselve s 



as computed from our astrometry with their cataloged values in UCAC3 (jZacharias et al.ll2009l ) 
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Given that the brightnesses and spectral types of the reference stars have approximately the same 
distribution for all our targets, which are also in the same general direction in the Galaxy, we can 
then find a correction of our CAPSCam proper motions to absolute proper motions. Leaving out 
TWA 13 and 15, which are visual binary stars whose proper motions in the UCAC3 are suspect, 
ten of our 14 stars have UCAC3 proper motions. CAPSCam-determined proper motions are indeed 
well correlated with the UCAC3 values, with a mean offset of -8.8 it 5.5 mas/yr in RA and 1.1 it 
5.1 mas/yr in DEC (Figured]). We correct these biases to find the proper motions of TWA sources 
without UCAC3 measurements, i.e. TWA 13AB, 15AB, 26, and 29. 



3. Results 

The results of our parallax survey are presented in Table [21 including the relative proper 
motions and parallaxes from the iterative solution, as well as the estimates of zero-point parallax 
in each field and the resulting absolute parallaxes. 

Two targets, TWA 13 and TWA 15, are visual binaries for which we obtained independent 
astrometry for the two stars in each system. Their solutions agree within their uncertainties. 



3.1. Notes on Individual Sources 

3.1.1. TWA 5 



TWA 5 (i.e. TWA 5A) has a companion brown dwarf (BD) TWA 5B (jWebb et al.lll999l ). It is 
visible in CAPSCam images taken during good seeing (Figure [2]), but the BD is faint enough and 
widely separated enough that it does not contribute to the PSF centroiding of the primary star. 

We measure the location of 5B in the 2009 April 11 and 2009 June 9 epochs, which both had 
excellent seeing of ~0'.'7. Because we were integrating on TWA 5A in the GW at 0.2 s and the 
reference FF for 30 s or 12 s, respectively, we have 2250 (720) GW images to use for "lucky imaging" 
to select the best few hundred images on each date. For each epoch, we shift and add these on 
the brightest pixel to form a final image with high quality - FWHM of 0'.'54 on 11 April and 0'.'60 
on 9 June. In these images, we measure the separation and position angles of the brown dwarf 5B 
with respect to 5A. The location of 5A is well-determined by the shift-and-add process. To find 
the TWA 5B centroid, we examine individually the flux in slices of each row and column around its 
location. In each slice, we subtract a smooth continuum from the bright star, and then fit the peak 
produced by 5B with a Gaussian. We then average the individual slice locations weighted by the 
height of the Gaussian in each slice to produce separate centroids in RA and DEC. Unfortunately, 
TWA 5B falls nearly on top of a diffraction spike that limits our ability to accurately centroid. The 
average separation on the two dates is 2'.'00 it OdO. 

We use 2MASS sources in the larger astrometric frame to solve for any PA offset between 
CAPSCam and the 2MASS reference frame and find offsets of 0.17 and 0.22 deg on the two dates. 
A match of the CAPSCam and 2MASS fields is also used to measure the CAPSCam pixel scale 
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Fig. 1. — Comparis on of proper motions of the ten TWA stars both reported in UCAC3 
([Zacharias et al.l |2009|) and measured in our survey. The dashed lines show the mean offsets of 
8.8 mas/yr in RA and -1.1 mas/yr in DEC between our relatve proper motions and the absolute 
UCAC3 ones. 
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of 0.1958 ± 0.0005 arcsec pixel"^. The average position angle of B with respect to A in the 2009 
epochs is 358.3 ifc 2.0 degrees. Both th e separation and position angle are consistent with the 
measurements of lNeuhauser et al.l (|20ld ) taken one year earlier. 



TWA 5A itself is a binary whose orbit of period 6 yr was measured with AO imaging (jKonopacky et al 



20071 ). We determine its distance to be 50.1 it 1.8 pc. The total dynamical mass of the components 
from the astrometric orbit depends strongly on distance (ocD^ ), which was tak en in Konopacky 
et al. to be 44 lb 4 pc from the convergent point analysis bv iMamajekl (j2005l ). Our parallactic 
distance is a 14% change from the Mamajek estimate, which makes a significant change in the 
inferred dynamical mass. The mass is now estimated at 1.05 ± 0.22 Mq instead of 0.71 ± 0.14 Mq. 

When we c orrect the absolute magn itude of TWA 5 A for binarity using Amag~l.l at H-band 



as measured by iKgnooackv et a 
agreement with the 



Baraffe et al 



(120071 ). we find that this brings the dynamical mass into good 



( 19981 ) pre- main sequence tracks, which predict a total mass of 1.2 
Mq and an age of 10 it 1 Myr from our interpolation in Section [4. II (see also Figure 4 of Konopacky 
et al.). 

This assumed that there were two and only two components in TWA 5A. However, there 
have been suggestions of additional rapid, high velocity variability that would suggest a separate 
spectroscopic binary within TWA 5Aa or 5Ab. ([Torres et al.ll2003l ) report that TWA 5A may be 
double- lined with possible day-scale velocity changes, but does not present individual RV measure- 



o 
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Fig. 2. — Image of TWA 5 A and 5B from the sum of 782 individual 0.2 s images taken 2009 April 
11 that were shifted to stack the brightest pixel at the center. TWA 5B is visible 2"north of the 
primary. 
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ments. iReidI (|2003l ) find two components separated by 30 km s^ , and may have observed the 
AO binary near periast ron. Periast rons were in 1998 .40, 2004.34, and 2010.28 according to the 
(jKonopackv et al.l 120071 ) ephemeris. ISong et al.l (j2002l ) report a single velocity of -30.6 it 6.6 km 
s~^ measured in 2001-2002 at apastron. Other measurements find mean velocity of ^13 km s 
with variabihty and line shapes consistent with motion of a few km s~^ ( Torres et al.l 120061 . 



2008 



Shkolnik et al. 



2011 



). TWA Aa and Ab are o f very similar spec t ral type an d brightness, and have 
significant rotational velocity > 30 km s~^ ( Torres et al.ll2003l : iReidI l2003l ) . so their lines should 
be difficult to distinguish. The AO orbit predicts a maximum relative velocit y of ^18 km s~^. I f 
there is a third star within TWA 5A, then the age of TWA 5Aa/Ab from the (IBarafFe et al.lll998l ) 
tracks would be underestimated (i.e. they would have lower luminosity). Their masses from the 
tracks would be essentially unchanged because the tracks are nearly vertical (see Figure [5]) , but 
the addition of the mass of the third component would imply that the tracks overestimate the 
dynamical masses of Aa/Ab alone. 



3.1.2. TWA lie 



Kastner et al.l ([200a) identified 2M1235 as a likely tertiary companion to HR 4796A, ie. as 



TWA lie. Its parallactic distance of 69.0 ifc 2.4 pc iii deed agrees well with the Hipparcos distance 
to HR 4796A of 72.8 pc ± 1.7 pc dvan LeeuwenlboO?! ^. 



3.1.3. TWA 16 



Zuckerman et al.l (|200ll ) reported that TWA 16 was a close visual binary with separation ~0'.'67 
and flux ratio ~0.9. The CAPSCam images reveal an elongated source in all epochs and resolve 
two sources during the observations with the best seeing (Figure [3|). The pipeline does not identify 
two sources there, however, so the measured parallax and proper motion are for the photocenter 
of the system. For the 2009 June 8 data, which had the best seeing, we used "lucky imaging" to 
select the best few hundred GW images. We shift and added these on the brightest pixel to form a 
final image with high quality and then used PSF fitting to measure the separation as 0'.'61 it 0'.'02, 
PA as 44°, and flux ratio as 0.97. 



Discussion 



4.1. Pre-Main Sequence Track Ages 



We determine absolute magnitudes using our parallaxes and 2MASS apparent magnitudes 



We a lso correct for binarity for six sources. TWA 2 is a visual binary with Amag~l (jWebb et al 
999l). TWA 5 A (see Section I3.1.ip is a close (speckle/ AO resolved) bin ary with Amaggil.l 



Konopacky et al.l 120071 ). TWA 14 is an approximately equal brightness SB (jjayawardhana et al . 
2006|) and Shkolnik, 2011, person al communication). TWA 16 is a visual binary with a flux ra- 

TWA 20 is a SB for which we assume the components are 



tio of 0.9 (Zuckerman et al 



2001 



_9_ 



equal brightness (iJayawardhana et al.ll2006l ). TWA 23 is a SB with equal brightness components 



([Shkolnik et al.ll201ll ). No attempt has been made to correct for extinction, which is small in the 
near-infrared due to the closeness of the stars and the absence of edge-on optically thick circumstel- 
lar disks in our sample. The tabulated uncertainties in absolute magnitude include the photometric 
uncertainty in the 2MASS H-band measurement and the parallax uncertainty. 

We obtain effective temperatures b y converting literature spectral types to ter nperature using 



the iii termediate scale of lLuhmanI (|l999l ) for the M-type stars and tabulated values in lHartigan et al 



(j 19941 ) for the earlier type stars. These are all given in Table [3l Most of the spectral types come 
from recent compilations that use TiO band or other spectral index fitting and should be mutually 
consistent and good to ~75 K. TWA 25 has no published spectral type; to obtain its temperature, 
we fit a Kurucz model to its photometry. TWA 29 (DEN1245) is the latest spectral type object 
for which we measured a parallax and sits near the M-L transition. The spectral-type to effective 
temperature conversion is not well known for such objects, and we approximate it at 2250 K. 
Finall y, we must iiote that historical optical an d new infrared spectral types for TW Hya do not 



(|Webb et al.lll999|;_ IVacca fc Sandellll201ll '). TW Hya's optical spectrum has been typed as 



Vacca &: Sandell 



(|201ll ) determination of 3400 K is likely to be too cool (N. 



agree 

K7V, 4000 K, but the 

Calvet, 2012, personal communication); we have chosen an intermediate value of the Tgff of 3615 
K. 

Theoretical isochrones overplotted with the data for all stars with parall axes (literature as well 
as this work) are shown in Figure HI We have chosen the lBaraffe et al.l (| 19981 ) tracks with Y=0.775, 
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Fig. 3. — Image of the close binary TWA 16 on 2009 June 8 with a separation of 0'.'61 it 0'.'02 and 
PA of 44°.' 
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mixing length parameter 1 for m<O.6M0 and Y=0.282, mixi ng length 1.9 for m >0.6Mf7^, based on 
their relative success in reproducing multiple star coevality (jWhite et al.lll999l ). The combination 
of the isochrones with different helium abundances creates a small temperature discontinuity at 
3500-3700 K, which is unfortunately the temperature range of many of our stars (Table [3l) , but we 
interpolate over this region anyway. We use the DUSTY isochrones (IChabrier et al.ll2000l ) for the 
stars with T^s <2900 K. We interpolate the combined theoretical tracks to estimate the ages of all 
the stars at their nominal positions in absolute magnitude-Tefr space. The median age is 10.1 Myr. 
Individual ages are also given in Table El along with all the data plotted on the tracks. Because 
of the abundance of stars with Teff~3600 K, Figure [5] shows a zoomed view of this part of the 
diagram. From the parallax uncertainty alone, the typical age uncertaint y on each star is 3 Myr . 
An age of 10 Myr is consistent with previous estimates, as summarized in [Fernandez et al.l (j2008l ). 



4.2. Galactic Space Motion and Membership 



For the analysis of the kinematics of the group, we use all stars with known proper motions 
(from UCAC3 if available), parallaxes, and radial velocities. These are given along with calculated 
UVWs in Table IH Uncertainties in all three quantities are p ropagated into the fi nal velocity 
uncertainties. We also use positions from the 2MASS catalog (ISkrutskie et al.ll2006l ). We then 
compute the mean velocity of the entire TW Hya association. 

The uncertainty weighted average velocities in Tableware [-10.1, -17.9, -8.0] it [0.2, 0.2, 0.1] 
km s~^. Two stars have velocities more than 3a from the mean of the association in at least 
two direct i ons: TWA 9 A and TWA 22. T WA 22 was already suspected not to be a member by 
Mama j eld (120051 ) and iTeixeira et al.l (120091 ). based on similar velocity arguments. However, TWA 
9A is a "classical" member used to define the convergent point in the Mamajek analysis. Given 
its discrepant age in Table [3l as well as its discrepant velocity, we conclude it is not a member or 
that its Hipparcos distance is underestimated. The new average velocities after excluding these two 
stars are: [-10.9, -18.2, -5.3] ± [0.2, 0.2, 0.2] km s'^ The standard deviation of the total velocities 
is 2.0 km s~^, and the RMS deviation in the total velocity from the mean total velocity is 1.9 
km s ~^. These are consider ably lower dispersions than obtained when photometric distances are 
used (jFernandez et al.ll2008l ). 



The uncertainties on the mean velocities above are computed assuming the stellar velocity 
distribution is Gaussian, i.e. standard deviation of velocities divided by the square-root of the 
number of stars. However, a K-S test reveals that the velocities are not Gaussianly distributed. 
Nor are they uniformly distributed between their minimum and maximum values. 

To trace the stars back in time, we take their present positions and three-dimensional space 
velocities and compute their locations in Galactic coordinates for time-steps back every 100,000 
years. To treat the distance and velocity uncertainties properly, we do this in a Monte Carlo for 
10000 trials selecting each star's distance and velocity in each direction randomly in each trial but 
distributed assuming the uncertainties for each individual star are Gaussian. Then, the centroid 3D 
location of the stars at each time is computed as well as the average distance of the stars from this 
centroid. The time of best convergence is defined to be when the average distance is minimized. We 
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Fig. 4. — TWA stars with parallaxes, and there fore a bsolute magnitudes, p lotted on theoretical pre- 
main sequence tracks from iBarafFe et al] ( 1998 ) and Chabrier et al. (200C). Effective temperatures 
largely come from converting published spectral types to temperature and have uncertainties (as 
shown) of at ~75 K (one half of a spectral type). The apparent scatter in ages is discussed in 
section 14.11 and Fig. El 
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Fig. 5. — Closeup of Figure H] in the region around Tcg=3600 K where most of the stars sit. 
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tested that the Galactic potential does not significantly affect the motions over the short timescale 
of 10 Myr. 

The present average distance from their mean location for the stars in Table [H excluding 9A 
and 22, is 20.4 pc. The present-day locations of the TWA stars with measured parallaxes and 
radial velocities is shown in Figure [H The velocities are nearly parallel, as shown in Figure El The 
nominal closest approach of all the stars is 2 Myr ago, but has a mean distance of 19.2 pc from the 
center, meaning that there is no time in the past when the stars are significantly more concentrated 
than they are today. At the mean age of 10 Myr established in section HTTj the mean distance from 
the center is 34 pc. 

Figure [8] displays a histogram of the ages excluding DAB, 22, and 29. There is a tail of stars 
to apparently larger ages while there are no stars with inferred ages less than 3 Myr. The best 
Gaussian fit to the age distribution has a mean of 9.5 Myr and standard deviation of 5.7 Myr. To 
assess the reliability of these values for such a small sample, we repeated the histogram fit in a 
Monte Carlo. In each trial, the age of each star was drawn from a Gaussian distribution based on 
that individual's star mean age and absolute magnitude uncertainty on the iBaraffe et al.l (|l998l ) 
tracks. The typical age uncertainty on each star is 3 Myr. The mean age over all the trials was 8.7 
Myr and the width of the age distribution over all the trials was 6.5 Myr. Thus, the age histogram 
is robust against the individual age uncertainties, and the width of the distribution compared to 
the typical age uncertainty indicates that there is a real spread in derived ages. 



14 



10 - 



Q 







lie 

16 
• 


' L 

15B_ 

15A 




*23 25 

/ 

28 


20 




9A 
2(? , 

13A 

4 
• 

2 


14 






21 

- 22 
I ■ 












Q 

N 



10 20 30 40 50 

X Distance (toward GC) [pc] 



60 



-20 



-60 



-100 



22 










• 












26 








4 


• 9A 








; 


5 ■ 25 








2 


21* . ^ 








% 


• • 23 28 








1,27 


l3/m 


• 








13B • 


lie • 


16 29 






12 


llA 


•. • 










20 
















14 
















15A ^ _ 












15B 

- 



10 20 30 40 50 

X Distance (toward GC) [pc] 



60 



Fig. 6. — Location of the TWA stars currently in Galactic coordinates. The size of the symbols 
is inversely proportional to age; that is, larger symbols indicate younger stars, as determined in 
Section [4. 1[ There is no apparent correlation between age and location. 
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5. Conclusions 



The identification of TWA members has largely been based on youth plus similarity to the 
young star TW Hya in terms of location on the sky, proper motion, and radial velocity. Co-evality 
was thought to follow under the assumption that these stars with similar motions formed from the 
same raw material at the same time. With parallaxes to fourteen primary stars identified as TWA 
members, we have greatly expanded the knowledge of the kinematics of these young stars. We find 
that although they do share a common space motion, the stars do not appear to have formed in a 
concentrated volume with a well-defined expansion velocity. The TWA stars appear to have formed 
over a larger volume than they presently occupy. 

Nor do the stars appear to be completely co-eval, as the stars studied here have ages that range 
from 3~23 Myr as derived from their locations on pre-main sequence tracks. Although uncertainties 
in the distance and the effective temperatures allow for several Myr of uncertainty in individual 
ages, it would be extremely difficult to force them all to a common age. This apparent age spread 
could be due to a real difference in the times of formation of the stars or it could be due to the 



lasting effects of episodic accretion (jBaraffe et al.ll2009l ) 



The spatial distribution of these nominal TWA stars of 40-60 pc is largely filamentary in nature, 
which naturally leads to some conclusions about their provenance. Perhaps the stars formed in an 
extended wisp of molecular cloud, probably one related to the Scorpius-Centaurus complex that 
is nearby on the sky but 80 pc further away. The Galactic V and W velocities of TWA are very 
similar to t he older Upper C entaurus Lupus (UCL) and Lower Centauru s Crux (LCC) subgroups 
of Sco-Cen (jchen et ahlboilh at ages of 15-17 Myr (jMamajek et ahlbood ). TWA is located near to 
LCC in Galactic coordinates but their separation in distance, age, and velocity distinguish the two 




Fig. 7. — Three-dimensional space motion of the TWA stars with measured parallaxes and radial 
velocities in Galactic coordinates. Black dots show the present locations of the stars and the colored 
points show the motion in 200,000 yr timesteps for 15 Myr. The parallel velocities mean that the 
stars are never much closer together than they are at present. This plot excludes TWA 9A and 22. 
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groups. iFernandez et alj (j2008l ) showed that TWA was sim45 pc from LCC 8 Myr ago and could 
have been subjected to 0.5 supernovae per Myr. Multiple supernova shocks could have triggered 
star formation in dense parts of the filamentary progeni tor TWA cloud ove r a few million years and 
not in a regular progression from one side to the other. (jOrtega et al.ll2009l ) also suggest that stellar 
winds and supernovae from LCC and UCL could compress gas in the region of TWA, although the 
mechanism for an extended time of star formation is less clear in this case. 



The Las Campanas Observatory staff and operators of the duPont telescope, particularly Oscar 
Duhalde, Javier Fuentes, Herman Oliveras, Patricio Pinto, and Andres Rivera, made the CAPSCam 
observations smooth and efficient. Rebecca Rattray made helpful analyses of portions of these TWA 
data during an undergraduate internship at DTM in 2009. CAPSCam was built with support from 
the NSF ATI program and Carnegie Institution of Washington. We acknowledge support for 
the observing by the NASA Astrobiology Institute under cooperative agreement NNA09DA81A. 
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Fig. 8. — Age histogram of the TWA stars with parallaxes, excluding DAB and 22, which have 
kinematics inconsistent with membership, and 29, which has a highly uncertain Tefr and therefore 
uncertain age. The bin size is 3 Myr, which is the median uncertainty in the individual ages based 
on the uncertainties in their absolute magnitudes (i.e., parallaxes). The best fit Gaussian to the age 
distribution is overplotted with the dashed line and has a mean of 9 Myr and standard deviation 
of 6.8 Myr. The width of the age distribution exceeds what would be expected for a population of 
stars formed in a single burst. 
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Table 2. Astrometric Results 



Target 




IJ'RAcosDEC,rel 


IJ'DEC,rel 


Zero-point 
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Table 3. Stellar Parameters 



Star Sp. Type Sp. Type Ref Teff H(abs) H unc BCAH98 age 



1 
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0.26 
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3 
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Table 4. Velocities 
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1.7 


-11.4 


3.8 


-19.8 


2.6 


-4.1 


1.9 


16 


12.8 


0.5 


-49.2 


-21.2 


1.6 


0.8 


9.0 


0.4 


-9.7 


0.8 


-18.6 


0.6 


-6.0 


0.4 


20 


12.9 


0.6 


-64.4 


-28.6 


3.1 


1.0 


8.1 ^ 


4.0 


-14.2 


2.3 


-21.0 


3.4 


-9.5 


1.3 


21 


18.2 


0.5 


-61.9 


15.0 


1.2 


1.4 


17.5 ' 


0.8 


-12.0 


0.5 


-20.2 


0.8 


-4.9 


0.4 


22 


57.0*^ 


0.7 


-175.8 


-21.3 


0.8 


0.8 


14.8 ^ 


2.1 


-8.0 


0.4 


-17.1 


2.0 


-9.0 


0.1 


23 


18.6 


0.5 


-72.7 


-29.3 


0.9 


0.9 


8.5 ° 


1.2 


-10.6 


0.6 


-18.2 


1.0 


-5.4 


0.6 


25 


18.5 


1.2 


-74.0 


-27.7 


0.8 


0.8 


9.2 ' 


2.1 


-10.7 


1.4 


-18.7 


1.9 


-5.6 


1.0 


26 


23.8 


2.6 


-89.9 


-26.5 


4.2 


2.6 


11.6 j 


2.0 


-10.7 


2.0 


-18.8 


1.9 


-3.8 


1.3 


27 


19. 0'^ 


0.4 


-62.7 


-22.8 


1.7 


2.8 


11.2 j 


2.0 


-9.1 


0.6 


-16.6 


1.8 


-6.7 


1.0 



Note. — Parallaxes are from this work unless otherwise noted in footnotes. Sources of RVs are given in footnotes. 
''Hipparcos - Ivan LeeuwenI ll2007h 
iTeixeira et al.l ll2008h 

'^Weighted average of lBiller fc Close! l l2007h . ICizis" et al. I ||2007|) . and lDucourant et al. I ||2008|') 
^Torres et al] ll2003h 
HShkolnik etal] ll201ll ') 

^Bright Star Catalog V/50 



f5 Assume same RV as for HR 4796B from lStauffer et al.l lll995h 
iReidI i200i) 
ISone et alj booA ) 
jMohantv et alj 1 I2OO3II 



